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The  NSWC  Aeroprediction  Code  has  been  extensively  Applied  to 
the  prediction  of  static  and  dynamic  aerodynamics  of  missile  con¬ 
figurations.  Major  extensions  have  recently  been  inade  to/the  code 
to  extend  Its  capability  to  07 £  M  s_< .8  and  Of  £(${<.  180  and  also 
to  improve  the  transonic  inviscid  "body  alone/ static^aerodynamic 
predictions  and  the  dynamic  derivative  predictions  for  all  Mach 
numbers.  The  theoretical  basis  for  the  code  extensions  are  out¬ 
lined  and  previous  methods  are  briefly  reviewed.  The  code  is 
evaluated  through  comparisons  of  computational  examples  with 
experiment  for  body  alone,  body-tail  and  body-tail-canard  con¬ 
figurations.  The  speed  and  accuracy  of  the  code  are  ideal  for  use 
in  preliminary  design.  Examples  of  design  applications  to  specific 
tactical  weapon  configurations  are  presented. 
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INTRODUCTION 


A  continuous  need  exists  for  estimating  the  aerodynamic  characteristics 
of  a  wide  variety  of  tactical  missile  and  projectile  configurations,  especially 
in  the  preliminary  or  conceptual  design  phase.  To  meet  this  need,  the  Navy 
(in  cooperation  with  the  Army)  undertook  the  development  of  a  rapid,  inexpen¬ 
sive,  easy  to  use  Aerodynamic  Prediction  Code  in  1971.  The  code  was  developed 
so  as  to  handle  fairly  general  wing-body-tail  configurations  and  hence  have 
direct  application  to  a  high  percentage  of  tactical  weapon  designs.  Preliminary 
versions  of  the  code  were  published  in  1972,  1975,  and  1977.  The  changing 
mission  requirements  for  both  current  and  future  weapons  has  dictated,  however, 
the  need  to  revise  and  extend  the  capabilities  of  the  1977  version  of  the 
Aeroprediction  Code,  which  was  limited  to  M ^  £  3.0  and  small  angles-of-attack 
£  15°),  to  higher  Mach  numbers  and  angles^-of-attack. 

The  objective  of  the  current  effort,  which  is  nearing  completion,  is  to 
extend  the  1977  version  of  the  Aeroprediction  Code  to  M  ■  8  and  a  »  180°. 

In  addition,  modification  of  some  of  the  existing  methods  due  to  advances  in 
the  state-of-the-art  and  computer  program  optimization  is  desirable. 
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The  general  approach  of  the  code  development  has  been  to  combine  existing 
and  newly  developed  computational  methods  into  a  single  computer  program.  The 
basic  method  is  that  of  component  superposition  where  the  body-alone,  lifting- 
surface-alone  and  interference  contributions  are  added  to  obtain  total 
configuration  aerodynamics.  The  code  development  has  occurred  in  four  increments. 
The  first  throe  of  these  increments  were  previously  reported,  and  led  to  the 
development  of  a  code  capable  of  determining  the  aerodynamic  coefficients 
for  axisymmetric,  non  air-breathing  configurations  with  up  to  two  sets  of 
lifting  surfaces  for  low  angles-of-attack  and  Mach  numbers  to  3.0.  The  results  of 
the  fourth  increment,  required  to  meet  the  stated  objective,  is  the  subject 
of  this  paper.  Program  plans  for  this  effort  and  so^e  early  results  were 
presented  at  the  11th  Naval  Symposium  on  Ballistics.  The  theo^igs  used, 
outlined  briefly  here,  are  discussed  in  more  detail  elsewhere.  *  The 
resulting  code  has  computational  times,  required  for  the  estimate  of  static 
and  dynamic  aerodynamic  coefficients  for  a  body-tail-canard  configuration  for 
one  freestream  condition,  that  are  in  CPU  seconds  on  a  CDC  6700  computer  as 
opposed  to  minutes  or  hours  often  required  for  more  detailed  physical  and 
numerical  models.  The  accuracy  obtained,  however,  is  compatible  with  that 
required  for  preliminary  or  intermediate  design  estimates. 

Numerous  computations  for  a  variety  of  configurations  have  been  attempted 
on  the  Extended  Aeroprediction  Code  and  the  results  compared  with  available 
data.  Representative  comparisons  with  experiment  used  in  evaluating  the  code 
and  sample  applications  illustrating  the  use  of  the  code  to  achieve  improved 
performance  will  be  discussed  in  a  later  section. 


CONFIGURATIONAL  GEOMETRY  AND  FREESTREAM  CONDITIONS 

The  most  complex  configuration  considered  is  illustrated  in  Figure  1. 

The  body  may  be  pointed,  spherically  blunted  or  truncated.  The  remainder 
of  the  body  may  consist  of  one  or  two  piecewise  continuous  nose  sections,  a 
constant  crossection  afterbody,  and  a  boattail  or  flare.  The  wing  or 
canard  fins  have  a  trapezoidal  planform  with  a  biconvex  or  modified  double¬ 
wedge  crossection  and  sharp  or  spherically  blunted  leading  and  trailing  edges. 
Tip  edges  are  assumed  parallel  to  the  freestream  at  zero  angle-of-attack. 

Fin  crossections  are  piecewise  similar  with  span.  No  camber,  twist,  dihedral, 
or  airfoil  distortion  is  considered.  Lifting  surface  sets  are  planar  or 
cruciform.  Horizontal, all-movable  control  deflections  in  the  plus  position 
are  considered.  Canard/wing  and  tail  fin  sets  are  aligned. 

For  various  Mach  number  and  angle-of-attack  regions  there  are  geometric 
restrictions.  These  will  be  elaborated  upon  in  later  sections. 

Freestream  condition  description  consists  of  Mach  number,  Reynolds 
number  per  foot  per  Mach  number,  and  angle-of-attack.  Roll  orientation  is 
considered  at  higher  angles-of-attack  only.  Inlet  and  exhaust  plume  effects 
are  not  considered. 


SUMMARY  OF  ANALYSIS  METHODS 

The  analysis  methods  will  be  discussed  in  general  terms.  For  more  detailed  ' 
theoretical  discussions  and  derivations  of  the  individual  methods,  the  reader  is 
referred  to  reports  currently  published  as  well  as  those  in  publication  (References 
2-7,  9). 
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BODY-ALONE  STATIC  METHODS 
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The  body-alone  low  angle-of-attack  static  prediction  methods  are 
summarized  in  Figure  2.  The  transonic  nose  pressure  drag  and  body  inviscid 
normal  force  methods  were  modified  under  contract  to  NEAR,  Inc.,  with  addition¬ 
al  modifications  accomplished  in-house.  The  second-order  shock-expansion 
extension  work  was  accomplished  under  contract  by  Professor  F.  DeJarnette  of 
North  Carolina  State  University. 

The  transonic/subsonic  nose  wave  drag  prediction  assumes  that  the  nose, 
boattail,  and  base  are  aerodynamically  isolated.  This  assumes  the  exis- 
tance  of  a  minimum  length  afterbody.  The  nose  shape  is  assumed  to  be  a 
spherically  blunted  tangent  ogive  or  a  spherically  blunted  cone.  To  obtain 
results  for  a  more  general  nose  shape,  linear  interpolations  between  the 
tangent  ogive  and  the  cone  pressure  drag  predictions  are  made  using  the 
initial  body  slope  at  the  sphere-nose  junction  and  the  nose-afterbody  slope. 
Thus,  for  a  zero  shoulder  slope,  the  tangent  ogive  value  is  computed  and  for 
equal  values  of  initial  and  final  slope  values,  the  cone  value  is  computed. 

The  tangent  ogive  pressure  drag  prediction  is  based  upon  interpolation  in 
a  table  of  values  of  Mach  number,  M^  ;  nose  length,  L^;  and  nose  spherical 
radius,  R^.  Ranges  of  values  are  .8  £  M^  £  1.2,  .75  <  Ljg  <  5.0  calibers, 
and  0.  <  Rjj  <  .5  calibers.  The  majority  of  the  pressure  drag  data  was  genera¬ 
ted  by  solving  the  Euler  equations  by  an  unsteady  implicit  time  asymptotic 
method.  A  portion  of  the  data  was  generated  by  solution  of  the  full  potential 
equation  by  the  method  of  South  and  Jameson.  Below  M^  =  .8  the  pressure  drag 
is  decayed  quadratically  to  zero  at  =  .5.  The  cone  pressure  drag  is 
obtained  from  a  blend  of  integrated  pressure  data  and  Taylor-Maccoll  solutions. 
The  pressure  drag  is  obtained  as  an  interpolation  in  a  table  of  values  of  M^ 
versus  cone  angle,  6  (bluntness  is  neglected).  Ranges  of  values  are  .5  < 

M  <1.2,  0  <  6  <  20°,  Below  M  =  .5  an  asymptotic  value  is  assumed.  The 
bo°attail  wave  drag  is  based  upon  small  disturbance  solutions  assuming  a  long 
afterbody  for  1.05  <  M^  <  1.2.  Below  M  ^  =  1.05  the  boattail  drag  is  decayed 
linearly  to  zero  at  M^  =  .95. 


The  transonic/subsonic  normal  force  prediction  is  for  0  <  M  <1.2.  The 
new  transonic  normal  force  prediction  is  based  on  a  combination  interpolation 
and  least  squares  curve  fit  in  Mach  number  and  geometric  parameters.  The  nose 
is  a  blunted  tangent  ogive.  Other  body  parameters  are  afterbody  length,  boat- 
tail  angle,  and  boattail  length.  Parameter  limits  are  0  <  Rm  <  .5  calibers, 

1.5  <  Lv|  <  5.0  calibers,  0  <  (afterbody  length)  <  5.0  calibers,  0  <  Lg 
(boattail  length)< 2.0  calibers,  and  0  <_  8g  <_  10°  for  the  conical  boattail 
angle.  Values  of  C^a  and  cMa  (about  the  nose)  are  obtained  by  solving  the 
Euler  equations  at  1°  by  an  unsteady  implicit  time  asymptotic  method  at  M^  = 
.75,  .90,  .95  and  1.2,  Currently  prediction  values  at  M^  =  .6  are  given  by 
the  earlier  model  and  at  M  *  1.2  by  a  low  supersonic  Mach  number  potential 
model.  Interpolation  provides  a  solution  for  .6  <  M^  <  1.2.  Below  M^  *  .6, 
the  earlier  model  was  used.  Currently  the  new  algorithm  is  an  input  choice 
to  the  program  since  neither  method  is  particularly  satisfactory  over  the  full 
range  of  conditions. 

The  small  disturbance  potential  solutions  for  the  low  supersonic  range  are 
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applicable  from  1.2  <  M^  <  M^J  M.,  a  program  input,  is  chosen  between  2  and  3.5 
depending  on  nose  shape  and  whether  or  not  lifting  surfaces  are  also  present. 


A  high  supersonic  prediction  method  is  used  for  M^  >  M^.  The  method  used  is  a 
modified  version  of  the  second-order  shock-expansion  method  which  predicts 
inviscid  static  coefficients  with  good  accuracy  for  bodies  with  short  after¬ 
bodies  and  a  flare.  The  wave  drag  predicted  is  adequate ^iotsnose-afterbody  ^ 
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configurations  but  is  poorer  for  bodies  with  a  boattail.  The  normal  force 
prediction  is  fair  for  nose-afterbody  configurations,  and  again  is  poorer 
for  bodies  with  negative  slopes.  The  pitching  moment  prediction  in  general 
is  poorer  not  only  for  the  higher  Mach  numbers  but  is  a  weakness  of  all  of  the 
methods  available  over  the  entire  Mach  number  range. 


LIFTING  SURFACE  AND  INTERFERENCE  STATIC  METHODS 

Lifting  surface  and  interference  methods  are  summarized  in  Figure  3.  The 
high  supersonic  methods  were  developed  by  Professor  F.  DeJarnette. 

Above  M  =  M  ,  a  shock-expansion  strip  estimate  is  used.  First  the  angle 
between  the  local^surfnce  normal  and  the  free  stream  velocity  is  found#  A 
local  oblique  shock  value  is  used  for  compression  angles  and  a  local  Prandtl- 
Meyer  value  is  used  for  expansion  angles.  The  pressure  distribution  is  indepen¬ 
dent  of  span  since  similarity  of  crossection  is  assumed.  For  a  blunt  leading 
edge  a  modified  Newtonian  distribution  is  assumed.  As  Mach  number  and  aspect 
ratio  increase  the  wing-alone  prediction  improves.  However,  interference 
effects  are  neglected. 

3 

The  methods  used  at  lower  Mach  numbers  were  previously  documented.  For 
the  symmetric  low  supersonic  drag  problem  no  swept  forward  trailing  edge  is 
permitted.  This  is  a  numerical  method  restriction  rather  than  a  physical 
restriction.  For  the  low  supersonic  lift  problem  no  subsonic  trailing  edge  is 
allowed.  In  addition,  the  zone  of  influence  must  not  include  the  opposite  tip 
edge  of  a  two  fin  planform.  For  the  drag  problem,  the  trailing  edge  sweep 
angle  is  cut  off  at  0°.  For  the  lifting  problem  the  Mach  number  is  kept  just 
above  or  at  the  critical  Mach  number. 


METHODS  FOR  COMPUTING  DYNAMIC  DERIVATIVES 

Methods  for  computing  dynamic  derivatives  are  summarized  in  Figure  4. 

For  the  fourth  increment  in  the  code  development  process,  improvements  and 
extensions  were  made  in  the  pitch  and  roll  damping  prediction  method. 

Dr.  L.  Ericsson  of  Lockheed  Missiles  and  Space  Company,  under  contract  to 
NSWC,  was  responsible  for  developing  a  method  to  improve  the  prediction  of 
pitch  damping  for  small  angles  of  attack.  In  the  low  supersonic  region  the 
method  is  restricted  by  the  slender  wing-body  analysis.  It  is  assumed  that  the 
slender-wing  body  analysis  is  usable  to  M  =2.  At  this  condition  the  aspect 
ratio  of  the  tail  is  restricted  to  less  tfian  2.3.  The  current  theory  is  res¬ 
tricted  to  body-tail  configurations.  For  the  body-alone  prediction,  when  the 
LMSC  model  deviates  too  greatly  from  the  older,  empirical  G.E.  SPINNER  code 
prediction,  the  SPINNER  prediction  is  selected.  Currently  the  use  of  the  LMSC 
algorithm  is  a  code  input  option.  For  body-canard-tail  configurations  only 
the  body-alone  prediction  is  utilized. 

For  the  earlier  methods  used  in  the  code  the  same  restrictions  on  Mach 
number,  which  prevent  trailing  edge  subsonic  conditions  or  the  opposite  side 
edge  lying  in  the  zone  of  influence,  which  were  made  for  the  lifting  surface- 
alone  normal  force  problem,  again  apply.  For  the  lifting  surface-body  problem 
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the  planform  considered  is  that  obtained  by  extending  the  leading  and  trailing 
edges  to  the  body  centerline. 

Improvements  were  also  made  in  the  method  of  computing  the  transonic 
empirical  roll  and  pitch  damping  for  configurations  with  lifting  surfaces.  At 
high  Mach  numbers  above  Mw  ■  M^,  a  strip  method,  based  on  the  static  normal 
force  strip  loading,  is  used  to  predict  roll  and  pitch  damping  for  configura¬ 
tions  with  lifting  surfaces.  Thus,  two  methods  are  available  for  predicting 
pitch  damping  for  the  entire  Mach  number  regime. 


HIGH  ANGLE-OF-ATTACK  STATIC  AERODYNAMICS 

The  method  used  in  the  code  is  a  direct  adaption  of  the  empirical  methodo¬ 
logy  of  Reference  7.  Restrictions  of  this  code  are: 

(1)  Mach  number:  .8  to  3.0. 

(2)  Angle-o f-attack:  0°  to  180°  for  isolated  components  (plus  position) 
and  0°  to  45°  for  body-tail  combinations  and  roll  angles  0°  to  180°. 

(3)  Tail.  Trapezoidal  plan  form,  edge  parallel  to  body  centerline 

(a)  Leading  edge  sweep  angle:  0  to  70  degrees 

(b)  Taper  ratio:  0  to  1. 

(c)  Aspect  ratio  (two  fins):  .5  to  2.0. 

(d)  No  control  deflection 

(e)  Tail  trailing  edge:  sweep  zero  and  parallel  to  body  base. 

(4)  Nose  length  (pointed  tangent  ogive):  1.5  to  3.5  calibers. 

(5)  Cylindrical  afterbody  6  to  18  calibers  long. 

(6)  Total  span  to  diameter  ratio  (two  fins):  1  to  3  1/3. 

Body  alone  and  body-tail  normal  force  and  Xcp  are  predicted  as  a  function 
of  Mach  number,  roll  angle,  and  angle  of  attack,  in  addition,  the  roll  moment 
is  also  predicted. 


ANALYSIS  EVALUATION  AND  SAMPLE  COMPUTATIONS  FOR  NEW  CODE  ELEMENTS 


This  section  will  present  examples  illustrating  both  the  strengths  and 
the  weaknesses  of  the  new  methods  adapted  and  integrated  into  the  NSWC  code. 
The  reader  should  remember  that  all  of  the  new  methods  met  our  accuracy 
requirements  in  general  over  the  range  of  applicability  although  in  some 
individual  cases  shown,  the  accuracy  was  poorer. 

An  illustration  of  the  capability  of  the  new  High  Mach  number  routine  is 
given  in  Figures  5-7.  Here  the  body  alone  static  aerodynamic  predictions  of 
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the  current  routine  are  compared  to  integrated  pressure  data  and  the  small 
disturbance  potential  codes^*®  for  a  blunted,  tangent  ogive  pose,  6.0 
caliber  afterbody  configuration.  Above  *  1.5,  the  comparison  of  the  High 
Mach  prediction  with  data  is  quite  good  for  the  forebody  drag,  C^,  (wave  drag 

plus  friction  drag)  while  the  normal  force,  C^a,  and  center-of-pressure,  X  , 

predictions  are  shown  to  degrade  somewhat.  In  general,  the  predictions  improve 
with  increasing  nose  length  and  Mach  number  (although  this  was  not  shown  for 

this  case  due  to  a  lack  of  data).  In  Figure  8,  comparisons  are  shown  for  a  NASA 

flared  body.  Small  modifications  were  necessary  to  the  body  geometry  to  make 
it  compatible  with  the  input  requirements  of  the  computer  code.  As  a  result, 
the  agreement  in  general  is  only  fair  except  for  C  at  the  higher  Mach  numbers 
where  the  agreement  is  good.  In  general,  the  C„  prediction  will  be  poorer  for 
bodies  with  negative  slopes  (i.e.  boattails) .  l?  is  also  shown  that  X 
is  rather  poorly  predicted  by  all  three  computational  methods.  C 

In  Figures  9  and  10  are  shown  comparisons  of  the  shock-expansion  strip 
theory  for  wings  with  the  low  Mach  number  small  disturbance  estimates.  For 
larger  aspect  ratios  and  higher  Mach  numbers  (smaller  zones  of  influence)  the 
comparison  is  improved.  No  interference  is  considered  for  M  >  M„ .  However, 

.  03  y  ' 

the  C  prediction  is  generally  on  the  high  side  which  partially  compensates 
for  this  omission  In  these  examples  A.  is  the  leading  edge  sweep,  A  is  the 
aspect  ratio  and  A  is  the  taper  ratio.  The  wing  crossection  is  thatwof  a  ' 
symmetric  diamond. 

Figures  11  and  12  illustrate  applications  of  the  code  to  body-canard-tail 
and  body-wing  configurations.  Data  for  the  TMX-1751  configuration  include  the 
contributions  of  body  strakes  (not  shown  in  the  sketch).  These  strakes  were 
not  included  in  the  model, due  to  limits  on  the  number  of  lifting  surfaces  and 
the  small  aspect  ratio .which  partially  accounts  for  the  C  and  X  differences. 
The  comparison  is  good  for  the  TMX-187  configuration  with  t:he  flarid  afterbody. 

Figures  13-16  illustrate  the  capabilities  of  the  adapted  Martin  High  angle- 
of-attack  empirical  algorithm.  The  slope  and  magnitude  of  the  C,,  predictions 
compare  quite  well.  However,  the  Xc  and  predictions  are  of  the  right  order 
of  magnitude  only,  except  at  fairlyC!igh  angles-of-attack  where  the  predictions 
improve  markedly  for  the  lower  Mach  numbers. 

As  previously  indicated,  the  modified  LMSC  dynamic  derivatives  model 
occasionally  experiences  a  total  breakaway  from  the  data  for  body  alone  con- 
figurations.  This  problem  is  illustrated  in  Figure  17  for  the  Army-Navy  Spinner 
configuration.  As  a  result,  the  current  combined  code  compares  the  LMSC  model 
predictions  with  the  G.E.  SPINNER  code  prediction.  If  the  deviation  is  large, 
as  in  this  case,  the  SPINNER  prediction  is  selected.  In  Figure  18,  the  relative 
capabilities  of  the  older  pitch  damping  prediction  method,  the  LMSC  model  and 
the  strip  theory  are  shown  for  the  Basic  Finner  configuration.  The  strip 
theory  is  seen  to  be  quite  adequate  at  the  higher  Mach  numbers  for  predicting 
both  the  pitch  damping  and  the  roll  damping  coefficient,  C  (see  Figure  19). 

The  overall  agreement  of  the  new  method  for  all  Mach  number!  is  considered  to 
be  good. 

The  relative  improvement  in  the  transonic  predictions  is  illustrated  in 
Figures  20-23,  In  Figure  20,  the  computed  transonic  nose  wave  drag  for  the 
M-117  Bomb  is  compared  with  data.  For  this  case,  the  NEAR  algorithm  is  shown 
to  improve  the  estimate  somewhat.  The  experimental  pressure  data,  however,  was 
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somewhat  sparse  for  an  accurate  numerical  integration.  The  improvement  is 
much  more  obvious  in  Figure  21  which  shows  a  comparison  of  the  NEAR  algorithm 
prediction  and  older  theory  with  data  for  blunted  tangent  ogive-cylinders 
with  6.0  caliber  afterbodies.  However,  there  is  no  apparent  improvement  in 
the  prediction  for  C„  and  X  in  the  transonic  flow  regime  as  shown  in 
Figures  22  and  23.  Neither  mlthod  seems  to  be  particularly  outstanding.  (A 
more  thorough  evaluation  is  present  in  Reference  5.) 


SPECIAL  CODE  FEATURES 


The  Aeroprediction  Code  which  applies  the  methodology  previously  indicated 
in  Figures  2-4  and  the  High  Angle-of-Attack  algorithm,  has  been  updated  to 
include  the  new  methodology  presented.  (Details  of  the  earlier  version  of  the 
code  are  available  in  Reference  9).  In  addition  to  the  new  computational 
methods,  other  elements  of  the  code  have  been  corrected  and  improved.  For 
example,  the  data  input  and  output  have  been  extensively  changed  to  simplify 
and  clarify  those  aspects  of  the  program  and  to  optimize  the  code  for  the  user. 
For  example,  the  wing  or  tail  input  data  consists  of  the  minimum  required  for 
the  computations  and  yet  allows  considerable  detail  to  be  given.  The  body- 
along  geometry  is  described  by  the  number  of  points  on  the  body  surface  for 
each  piecewise  continuous  segment,  a  single  logic  variable,  and  body  surface 
coordinates.  The  description  is  basically  independent  of  Mach  number  with 
only  minor  exceptions.  The  program  has  been  segmented  in  order  to  minimize 
the  storage  required  for  loading.  Currently,  the  code  requires  160K  octal 
central  memory  or  less.  For  more  universal  use,  all  FORTRAN  statements  will 
be  standard  ANSI  format  or  machine  independent.  For  missile  design,  the  output 
has  been  expanded  to  include  the  pitching  moment  coefficient  components  about 
the  given  moment  center  and  the  normal  force/pitching  moment  dependencies  on 
angle-of-attack  and  control  deflections.  Thus,  for  the  entire  configuration, 
the  normal  force  coefficient  for  a  given  Mach  number  can  be  described  as 
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where  6  and  6  are  the  canard  and  tail  control  deflections,  respectively,  and 
CNaa  t*ve  bo^  viscous  crossflow  term.  The  expression  is  similar  for  the 
moment  coefficient.  Utilizing  these  outputs  one  can  obtain  the  hinge  moments 
and  other  coefficients  needed  for  a  linear  aerodynamic  performance  analysis 
of  a  missile  configuration. 


APPLICATIONS  TO  DESIGN 

The  Aerodynamic  Prediction  Code  has  been  extensively  applied  to  the  design 
of  both  conventional  and  unconventional  tactical  weapons.  To  conserve  space, 
only  a  few  specific  examples  are  shown  here  which  hopefully  will  provide  some 
indication  of  the  code's  flexibility.  Many  more  examples  will  be  given  in  the 
design  manual  to  be  published  later  this  year  (Reference  6). 
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PROJECTILE  PERFORMANCE  OPTIMIZATION 

The  objective  of  this  study  was  to  optimize  the  performance  of  a  full- 
bore,  axisymmetric  projectile  to  obtain  the  minimum  time-of-f light  and  highest 
terminal  velocity  at  a  given  range  with  no  degradation  in  accuracy.  This  was 
to  be  accomplished  by  optimizing  the  shape  to  minimize  the  total  drag,  through 
determination  of  the  best  weight  to  obtain  the  optimum  ballistic  coefficient 
and  by  judicious  selection  of  the  center-of-gravity  location.  An  optimization 
scheme  developed  by  Hagar,  et  al.^^  (which  has  since  been  upgraded  by  Mogall) 
indicated  that  for  a  five  caliber  body,  the  minimum  drag  profile  would  be  simi¬ 
lar  to  that  shown  in  the  top  corner  of  Figure  24.  Methods  to  quickly  predict 
the  aerodynamics  of  this  profile  are  more  limited  than  one  might  at  first 
suppose,  particularly  with  respect  to  the  two-third  power  law  nose  contour  and 
the  1.8  caliber  boattail.  For  instance,  the  range  of  applicability  of  the 
DATCOM  is  for  boattail  lengths  greater  than  2.0  calibers  and  tangent  or  conical 
nose  shapes.  The  G.E.  SPINNER  code  includes  data  for  some  secant  ogive  nose 
shapes  which  can  closely  approximate  the  two-third  power  law  nose  but  is  limited 
to  boattail  lengths  of  less  than  1.0  calibers.  The  small  afterbody  length 
(approximately  .05  calibers)  also  restricted  the  use  of  additional  methods  and 
data.  The  Aeroprediction  Code  was  utilized  to  obtain  the  static  aerodynamic 
coefficients  for  this  configuration  so  initial  estimates  of  its  stability 
could  be  made.  Using  these  initial  estimates,  the  design  was  fabricated  and 
successfully  flown  in  the  30  MM  caliber  size.  Comparisons  of  the  theory  with 
the  ballistic  range  data  later  obtained  for  the  design  are  given  in  Figure  24. 
Good  agreement  was  obtained  in  each  case.  As  a  result,  the  predicted  values 
could  then  be  used  to  conduct  a  stability  analysis  throughout  the  trajectory 
with  greater  confidence.  Use  of  the  Aeroprediction  Code  therefore  allowed 
the  designer  to  eliminate  expensive  wind  tunnel  tests  and  the  bulk  of  ballistic 
range  tests,  even  though  extensive  changes  in  the  design  were  made,  and  demon¬ 
strate  the  benefits  of  the  proposed  design.  The  improvements  in  this  case 
were  substantial  as  shown  in  Table  I  below. 


Table  1.  Benefits  of  the  30  MM  Optimal  Projectile* 


30  MM 
Standard 

30  MM 
Optimal 

Improvement 

Range  (KM) 

3.0 

3.0 

— 

Average  CQ 

.36 

.22 

39% 

Time-of-Flight  (sec) 

5.70 

3.86 

32% 

Terminal  Velocity  (FPP) 

967 

1559 

61% 

Accuracy  (rad.  std.  dev.,  mts) 

.6 

.65 

— 

*  Computer  generated  trajectory  based  on  experimental  drag  data. 
STRUCTURAL  ANALYSIS 

An  example  of  application  toward  structural  design  was  presented  at  the 
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11th  Navy  Symposium  on  aeroballistics. 1  As  indicated  in  the  previous  paper,  the 
Aeroprediction  Code  was  modified  to  compute  the  total  sectional  normal  force 
coefficients  along  the  body  as  a  function  of  Mach  number  and  angle-of-attack. 

The  results  were  then  used  in  the  structural  analysis.  The  results  of  the 
analysis,  repeated  here  in  Figure  25,  indicate  that  structural  failure  of  the 
two  configurations  examined  could  occur  during  maneuvering  flight  due  to  the 
large  aerodynamic  bending  moment.  Flight  tests  confirmed  the  prediction  and 
the  design  was  successfully  modified. 


MISSILE  CONTROL  DESIGN 

The  use  of  the  code  in  preliminary  missile  design  was  recently  illustrated 
for  the  Advance  Point  Defense  System  (APODS)  missile  concept.  The  guided  APODS 
missile  is  currently  an  18.0  caliber  design  with  a  2.0  caliber  Van  KArmAn  ogive 
radome,  a  16.0  caliber  afterbody  and  four  cruciform  tail  fins  for  control  (see 
Figure  26).  In  order  to  obtain  an  estimate  of  its  performance  capabilities  and 
determine  the  optimum  control  gains,  the  Aeroprediction  Code  was  utilized  to 
provide  the  static  aerodynamics  at  angle-of-attack  and  the  moment  contributions 
relative  to  the  center-of-gravity.  The  estimated  static  and  dynamic  aerodynamic 
coefficients  for  the  APODS  missile  are  given  in  Figure  27.  These  results  were 
fed  into  the  computer  guidance  model.  A  sample  trajectory  and  the  tail  control 
autopilot  response  are  shown  in  Figures  28  and  29,  respectively.  Results  such  as 
these  can  then  guide  the  design  in  modifying  the  configuration  or  the  control 
gains  to  maximize  performance. 


SUMMARY 

An  Extended  Aerodynamic  Prediction  Code  for  rapid,  approximate  estimates 
of  the  static  and  dynamic  coefficients  for  guided  and  unguided  tactical  weapons 
has  been  developed.  The  range  of  applicability  of  the  code  is  0  Mb  £6.0 
(up  to  Ma  *  8.0  for  certain  body  alone  cases)  and  0°  £  a  £  45°  (up  to  a  *  180° 
for  computing  the  drag  on  certain  body  alonl  configurations).  A  large  number 
of  configurations  can  be  accurately  modeled  on  the  program.  Final  verification 
of  the  predictions  is  nearing  completion.  The  code,  accompanied  by  a  Design 
Manual  and  a  User’s  Guide,  is  scheduled  for  release  in  the  fourth  quarter  of 
FY81. 
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3.  New  Methods  for  Computing  Wing-Alone 
and  Interference  Aerodynamics 
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Figure  4.  New  Methods  for  Computing  Dynamic  Derivatives 


1-50 


O  DATA  (AFATL) 

-  AEROPREDICTION 

—  —  IMPROVED  Van  DYKE 
“ —  HIGH  MACH  CODE 


a-c r 


(a)  Lw  =  1.859,  Rm  --  0.375,  a  =  0 


(b)  Ln  =  2.357,  Rn  =  0  250,  a  =  0 


Figure  6.  Blunted  Tangent  Ogive-Cylinder, 
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Figure  8.  NASA  Flared-Body  Data  Comparison 
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Figure  11.  Comparison  of  Static  Aerodynamics  for 
Finned  Body  at  High  Mach  Numbers  (a  =  0°) 
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Figure  12.  Aerodynamics  for  TMX-187 
(M  =  4.65,  a  =  0°) 
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Figure  13.  CL.  Comparison  for  Air  Slew  Demonstrator 
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Figure  14.  X  Comparison  for  Air  Slew  Demonstrator 
Forward  of  Body  Midpoint 
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Figure  15.  C.  Comparison  for  Ait  Slew  Demonstrator 


Figure  17.  Army-Navy  Spinner  Pitch  Damping  Comparison 


Figure  18.  Basic  Finner  Pitch  Damping  Comparison 
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Figure  19.  C  Comparison  for  Army-Navy  Basic  Finner 
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Figure  20.  M-117  Bomb  Nose  Wave  Drag  Comparison 
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Figure  22.  CNa  and  Xcp  Comparison  for  Blunted  Tangent  Ogive-Cylinders 
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Figure  23.  Comparison  of  Theory  and  Test  Data  for  155  mm  Projectile 
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Figure  27.  APODS  Predicted  Aerodynamics  Using  Aeroprediction  Code  (a  -  15°) 


